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Abstract. We give a Riemann-Hilbert approach to the theory of matrix orthogonal poly- 
nomials. We will focus on the algebraic aspects of the problem, obtaining difference and 
differential relations satisfied by the corresponding orthogonal polynomials. We will show 
that in the matrix case there is some extra freedom that allows us to obtain a family of lad- 
der operators, some of them of 0-th order, something that is not possible in the scalar case. 
The combination of the ladder operators will lead to a family of second-order differential 
equations satisfied by the orthogonal polynomials, some of them of 0-th and first order, 
something also impossible in the scalar setting. This shows that the differential properties 
in the matrix case are much more complicated than in the scalar situation. We will study 
several examples given in the last years as well as others not considered so far. 
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1 Introduction 



The theory of matrix orthogonal polynomials on the real line (MOPRL) has its foundations 
in the seminal papers of Krein |421 143| (see also their account in the book of Berezans'kii [2]). 
For further historical background and analytic results, the reader is referred to the survey 
and to Chapter 4 of [51]. In many aspects the MOPRL resemble their scalar counterparts, 
especially where the proofs are based on basic properties of Hilbert spaces. Nevertheless, the 
non-commutativity of matrix multiplication and the existence of non-zero singular matrices add 
features to the theory that make MOPRL an interesting object of study. Moreover, many prob- 
lems for scalar polynomials are better understood or recast in terms of some matrix polynomials, 
see, for instance, [26] . 

Recall that a matrix polynomial of degree < n in C^^^ and a scalar variable x can be defined 
as an expression of the form 

Anx"- H h Aix + Ao, 
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where Aj^s are constant matrices in C^^^. In what follows we consider N fixed, and denote 
by P„ the family of matrix polynomials of degree < n in C^^^, as well as by P := IJn>o^"' 
We use preferably boldface letters to denote matrices, and standard font for scalars. We also 
use In for the N x N identity matrix, omitting the explicit reference to its dimension when it 
cannot lead anyone into confusion. 

In this paper we consider some aspects of the MOPRL theory, assuming that the orthogonality 
is given by an absolutely continuous measure on the line. More precisely, our starting point is 
a weight W = (Wij) : {a,b) — >■ GL{N,M), defined and positive definite on a finite or infinite 
interval (a, b) C M. We will assume that all Wij and W^j have finite moments: 

rb i-b 

I \x\^W{x)dx <oo, I \x\^W'{x)dx <oo, neNo:=NU{0}, 

J a J a 

where the integration of a matrix function is applied entry-wise. 

For any two -P, Q G P, the weight W induces two matrix- valued "inner products", 



and 



{P,Q)w= [ Pix)W{x)Q*{x)dx, 

{P,Q)w= [ Q*{x)W{x)P{x)dx = {Q\P*)w, 

J a 



where the asterisk denotes the conjugate transpose (or Hermitian conjugate) of a matrix. Due 
to this connection between both inner products, we restrict our attention to (•, We define 
also the norm 



\P\\w ■■= (Tr(P,P)w) 



1/2 



and assume that W is non-trivial, in the sense that ||P||t^ > for every non-zero matrix 
polynomial P. In this case (see jlH Lemma 2.3] or |5H Proposition 4.2.3]), {P,P)w is non- 
singular for every non-zero polynomial P, and we can easily implement a matrix analogue of 
the Gram-Schmidt orthogonalization procedure, which yields a unique sequence {Pn)n of monic 
orthogonal polynomials such that Pq = I, 

n-l 

P„(x) = x"J + ^ a^jx^ {Pn,Q)w = for every Q € Pn-i, n e N, (1.1) 

j=0 

as well as matrix polynomials {Pn)m "orthonormal" with respect to W, such that 

Pn{x) = KnPn, {Pn,Pm)w = Sn,ml- (1-2) 

Obviously, PnS are determined up to a unitary left factor, so we can speak about an equivalence 
class of orthonormal MOPRL, corresponding to the weight W. Additionally, if S is a constant 
non-singular matrix, then 

W{x) = BW{x)B* 

is also a weight, and (-Pn)n = {PnB~^)n is the corresponding sequence of orthonormal polyno- 
mials. Hence, more than a single weight we consider an equivalence class of weights given by 
its representative W. In particular, without loss of generality we can assume that W{xq) = I 
at a point xq € (a, 6). If this class contains a diagonal matrix- valued function, we say that W 
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reduces to scalar weights. A characterization of this fact, as shown in [23], is the commutativity 
condition 

W{x)W{y) = W{y)W{x), x,y £ [a,b]. 

This situation is considered trivial, and usually omitted from consideration. 

Work in the last few years has revealed a number of explicit families of MOPRL; in many 
cases they are joint eigenfunctions of some fixed differential operator with matrix coefficients 
independent of the degree n of the MOPRL. This study was initiated in [21], but nontrivial 
examples had to wait until [23^ [29l [34t [36] . A solution to the classification problem for those W 
whose MOPRL are common eigenfunctions of some fixed differential operator remains elusive. 
There are by now two methods that have yielded nontrivial examples: a) the connection with 
matrix- valued spherical functions for symmetric spaces |34[ [36] , and b) a combination of classical 
methods and some Lie algebra tools |23[ [24] • A necessary condition in terms of moments is 
given in [23] . and a necessary and sufficient condition in terms of the "ad-conditions" is given 
in [38]. In this approach one gets an explicit formula for one differential operator. Concerning 
applications, the recurrent relations have been used in the study of certain quasi-birth-and-death 
processes [301 ED 1321 [33] , in which case the differential operator plays no role. Recently, new 
applications of these processes have been related to urn or Young diagrams models [37] . 

A big boost in the research of the scalar OPRL has been their Riemann-Hilbert (RH) cha- 
racterization, introduced in the seminal paper of Fokas, Its and Kitaev [28], and complemented 
with a non-linear steepest descent analysis in a series of works of Deift, Zhou and collabora- 
tors [12[ [T3| [161 117j . This combination has allowed to establish extremely strong asymptotic 
results with applications in random matrix theory, in particular for unitary invariant ensem- 
bles [mill], determinantal point processes [8l[9l[lll[l5], orthogonal Laurent polynomials [36] [TT]. 
Painleve trascendents [I0l[4l], the Toda lattices [15], to mention a few. 

The RH characterization, even without the steepest descent analysis, allows one to prove other 
algebraic and analytic properties of orthogonal polynomials, as it was illustrated in [13] for OP 
on the unit circle, and in [401 Chapter 22], for classical families of OPRL. One of the goals of this 
paper is to extend these considerations to MOPRL, something that, surprisingly enough, has not 
been explored in depth so far. We will show that the RH formulation, which has a very natural 
block-wise generalization to the MOPRL case, allows one to reveal some subtleties hidden in 
the matrix case. For instance, we will derive identities and difference-differential equations that 
were unknown even for some of the explicit families mentioned above. The RH technique can 
be also a natural method for obtaining the differential relations for MOPRL with respect to 
general weight matrices. 

There have been several papers considering the block-wise RH problem of the flavor similar 
to the one discussed in this paper, but for singular (namely, rank 1) weight matrices. For these 
weights the MOPRL are connected with a class of multiple orthogonal polynomials that find 
applications in the analysis of determinantal point processes and non-intersecting stochastic 
paths (see e.g. the works of Kuijlaars and collaborators [8lll8 [ [T9 ] l45j) or of the multicomponent 
2D Toda lattice hierarchy [1]. 

This paper deals with strictly algebraic consequences of the RH formulation of matrix poly- 
nomials orthogonal with respect to a weight matrix supported on an interval [a, b] C M; the 
asymptotic analysis via the Deift-Zhou non-linear steepest descent method will not be discussed 
here. 

In Section[2]we will discuss two dual RHPs, which are uniquely solved in terms of the MOPRL. 
Standard arguments related to RHP allow us to obtain the three-term recurrence relation and 
the differential identities (ladder operators) for MOPRL. Both yield also the so-called Lax pair 
for the polynomials, which is an over-determined system (j2.28p . whose compatibility conditions 
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render further properties of MOPRL. All the results in this section have been previously obtained 
using different approaches, but the RH method will provide new and straightforward proofs. 

In Section [3] we will consider a special transformation of the original RHP, based on factori- 
zations of the weight matrix of the form W = TT*. The main goal is to obtain a RHP with 
constant jumps (independent of z), in order to simplify the differential relation obtained earlier. 
For simplicity, we will focus on the case where the support of the weight matrix is IR. The 
extension of the method to weights supported on semi-infinite or finite intervals is simple, see 
a short discussion at the end of this paper. 

We exploit also the non-uniqueness of the weight factorization, which in the matrix case 
gives us some extra freedom and yields non-trivial relations. In particular, a family of ladder 
operators, some of them of 0-th order, is obtained, a phenomenon that is not possible in the 
scalar situation. Furthermore, by combining appropriately the family of ladder operators, we 
will get the 0-th, first and second order differential equations satisfied by the MOPRL. Some 
examples of lower order differential operators appear in [5l [6} [35]. 

The considerations so far have been completely general, applicable to any family of MOPRL. 
In Section S] we narrow the analysis of Sections [2] and [3] to some relevant examples of MOPRL 
supported on M, obtaining new results even in the cases studied previously in the literature. 

In the final stage of preparation of this manuscript we learned about the preprint [1] which 
also uses the Riemann-Hilbert approach for the analysis of the matrix orthogonal polynomial. 
Although there is some overlapping between the results contained in our Sections [21 [3l and in [U 
Section 2], the focus of both contributions is different, and in this sense, complementary. 



2 The Riemann-Hilbert problem for MOPRL 

2.1 Formulation and basic properties 

In this section we discuss the Riemann-Hilbert problem (RHP) related to MOPRL with respect 
to a X weight matrix W supported on an interval [a,b] C M; this interval can be either 
bounded or unbounded. We assume W continuous and non- vanishing on (a, 6), and that at any 
finite endpoint of the support the weight W has at worse a power-type singularity, that is, W 
is of the form 

W(z) = \z-c\^-W{z), 7c>-l, (2.1) 

where c S {a,b}, c ^ =boo, and W is a bounded, continuous and non-vanishing at z = c. This 
class of weights comprises all the examples considered so far in the literature. 

As a convention, in what follows we write the 2N x 2A^ matrices partitioned into N x N 
blocks, as in p.2p below. Recall that A* stands for the Hermitian conjugate of the matrix A, 
as well as A~* = (A*)~^. We adopt the convention that for any matrix-valued function P(z) 
of a complex variable z, P* denotes the matrix-valued function obtained as 

P*(z) := {P{-z)r. 

The RHP for MOPRL with respect to a weight matrix W consists in finding a matrix function 
: C ^ (C2Afx27V g^jpjj ^^^^ 

(Yl) y is analytic in C \ [a, h\. 

(Y2) "K" has on (a, 6) continuous boundary values Y]^ (resp., Y^) from the upper (resp., lower) 
half plane, such that 



r"(x) = r"(x) 



In 




W{x) 
'"In" 



xe{a,b). (2.2) 



Properties of Matrix Orthogonal Polynomials via their Riemann-Hilbert Characterization 5 



(Y3) As z — )■ oo, for every m S N we have 



1=1 














(where the asymptotic term ©(l/z*""*"^) depends on n). 
(Y4) As z — )■ c, c G {a, 6}, c 7^ ±oo, we have 



Y^{z) 



Oil) 


0[h{z)) 


0{l) 


0{h{z)) 



where 



\z - c|'''^ if - 1 < 7c < 0, 
h{z) = { log \z-c\, if 7c = 0, 
1, if7c>0. 



(2.3) 



(2.4) 



Remark 2.1. Usually the asymptotic condition at infinity (|2.3p is stated for m = 1; however, 
it can be proved that (j2.3p for m = 1 implies that this condition holds for every m € N. 

Along with the RH problem (Y1)~(Y4) we can consider the following dual problem: finding 
a matrix function j/"' : C — )• C^^^^^ such that 



(yl) is analytic in C \ [o, h]. 

'In 


(y3) As z — >■ oo, for every m E N we have 



(y2) yl{ 



y"i{x) when x G {(i,b). 

In J 



y'^iz) 



z-'^In 

6 





"z^In 



(2.5) 



(y4) As z — >■ c, c G {a, 6}, c 7^ ±oo, we have 



y'^iz) 



\ 0(1) 
with h defined in (IZil). 



0{h{z)) 
Oil) 



It will turn out that (yl)-(y4) is related to the inverse of the solution of (Y1)-(Y4). 
For any integrable N x N matrix- valued function F on [a,b], 



C{F)iz) :- 



Fit) 



2Tri ./„ t — z 



dt 



defines the Cauchy or Stieltjes transform of F, which is a matrix-valued and analytic function in 
C \ [a, b]. Let us introduce the matrix polynomials of the second kind (of degree n — 1), defined 
by 



Qnix) 



Pnit) - Pnix) 



t — X 



Wit)dt, n > 0. 



(2.6) 
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We have 

2T:iC{PnW){z) = Qn{z) + 2mPn{z)C{W){z), 

where 2'KiC{W) is the m-function of the weight matrix W . 

Finally, if k„, is the leading coefficient of any corresponding normalized polynomial Pn, we 
denote 



(2.7) 



Observe that a priori 7„ depends on the selection of Kn- 

The following is a complete analogue of the well-known theorem from |28| , although its proof 
requires some additional considerations: 



Theorem 2.2. The unique solution of the RHP (Y1)-(Y4) is 



Y^iz) = Rn{z)Y\z) 



n > 0, 



where 



Y^{z) 



In 




C{W){z) 
In 



and the transfer matrix Rn is a matrix polynomial given by Rq^z) = I, 



Rn{z) 



^n'^(f) 

-27riK*_iP„_i(z) 



{2'ninn) ^Qn{z) 

— K.^„lQn-l(z) 



n G N. 



Analogously, the unique solution of the RHP (yl)-(y4) is 
= yO(z)r„(z), n > 0, 

where 



In 




-C{W){z) 
In 



and the transfer matrix Vn is a matrix polynomial given by ro(z) = /, 
rn{z) 



( -g;_i(z)Kn-i 

I 2-KiP*^-^{z)nn~i 



n € 



Moreover, for all n > 0, 

detr"(z) = 1, for all z £ C. 
Here denotes the transpose of the matrix A 






In 


-In 






{Y^{z))-\ 



(21 



(2.9) 



(2.10) 



(2.11) 



(2.12) 
(2.13) 



Remark 2.3. In the scalar case, when Y"' is a 2 x 2 matrix, there is no need to consider simul- 
taneously both RHPs (Y1)-(Y4) and (yl)-(y4) since the existence of {Y"-)'^ and ([ZTTD . ([232]) 
follow directly from (12.13p . In the general (2A^) x (2A^) case more care should be put in the 
analysis of the local behavior at the endpoints a, 6, and ()2.12p is no longer straightforward. 
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Remark 2.4. The solution of (Y1)-(Y4) satisfies the following symmetry relation: 
Viz) ■- 








i 


-In 



In\ 



(2.14) 



which yields some obvious consequences for the block entries of Y"". This relation is established 
using the invariance of (Y1)-(Y4) by such a conjugation. 



Remark 2.5. Alternatively, we can write the solution (|2.8p as 
Y^iz) = 



Pn{z) ^S^^^M 
-27ri-fn-lPn-l{z)\-27rijn-lC{Pn-lW){z) 



n eN, 



(2.15) 



where Pn denote the monic MOPRL of degree n for the weight matrix W. In the same vein, 



{Y^)-\z) 



-27TiC(WP:^,)iz)'rn-l\ 


-CiWP:){z) 


27riP„*_i(z)7„_i i 





n G N. 



(2.16) 



Remark 2.6. Although fv„ is defined up to a left unitary factor, the matrix coefficient 7„ 
in ()2.7p is unique. This is a consequence of the uniqueness of the solutions of the RHP above. 

Remark 2.7. Formulas (j2.8p and (j2.10p show that generically the behavior of Y'^ and (1^")~^ 
at the endpoints (and in general, any singular point) of the support of W is given by the local 
behavior of the Cauchy transform of the orthogonality weight, also known as its m- function, see 
[501 § 1-2]. 



Proof of Theorem The fact that ([22]) or ([2351) is a solution of the RHP (Y1)~(Y4) is 
established following the proof for the scalar case, see e.g. [12], and taking into account that the 
orthogonality of P„ in (jl.ip is equivalent to the homogeneous system 



J a 



x^Pn{x)W{x)dx = 0, 



i = 0, 1, . . . ,n - 1. 



The same applies to (fZTOD or ([ZTBD and the RHP (yl)-(y4). 

Consider the function Z^{z) = Y^{z)y^{z); from (Y2) and (y2) it follows that it has no 
jump across (a, 6), and by (Y4), (y4), 

Z"(z) = 0(/i(z)), z^ce{a,6}. 

Hence, has only removable singularities at the finite endpoints of the support of the weight, 
and thus is an entire function. It remains to observe that Z'^{oo) = I2N to conclude that 
Z^{z) = I2N for all z ^ C, which proves that y^{z) = {Y^{z))~^, as well as the uniqueness 
of both solutions. The first identity in (|2.12p can be established by direct calculation or by 
observing that this transformation carries the RHP (Y1)-(Y4) to the RHP (yl)-(y4). 
Finally, the scalar function detl^"'(z) is analytic across [a,b], and by (|2.8p . (|2.9p . 

detr"(z) = det Rn{z) 

can have only removable singularities at the (finite) endpoints of the support of W. Hence, 
detl^"(z) is an entire function; since by ()2.3p . detl^"(oo) = 1, we conclude that it is identi- 
cally 1. ■ 
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Proposition 2.8. Let Pn = i^nPn, n >0, be a sequence of orthonormal MOPRL, and let {Qn)n 
be the corresponding matrix polynomials of the second kind (|2.6|) . For these polynomials, define 



(2.17) 



Then 



Qn{z)P:{z)=Pn{z)Ql{z). 



(2.18) 
(2.19) 



Moreover, P*_i{x)AnPn {x) and Ql^^{x)AnQ n{x) are Hermitian for all x G M and n > 0. 



Remark 2.9. Identities ()2.18p and (j2.19p are also known as the Liouville-Ostrogradski formula 
and the Hermitian property, respectively. Both were originally derived for MOPRL in [20] , using 
a different approach. 

Proof. Using the explicit expressions for and (1^")~^ (written in the form (|2.8p and (|2.10p ). 
from block entries (1,1) and (2,2) of the identity Y^{Y^)~^ = I we obtain the Liouville- 
Ostrogradski formula (j2.18p while from block entries (1,2) and (2, 1) we get the so-called Her- 
mitian property (|2.19p . 

Block entries (1,2) and (2, 1) of (Y'^)~^Y^ = I yield the commutativity relations 

p:_,{z)AMz) = P:{z)A:Pn-i{z), Q:,„,{z)AMz) = Q:{z)AlQn_i{z), 

which show that P*_^{x)AnPnix) and Q^_i{x)AnQnix) are Hermitian for all x G M and 
n > 0. ■ 

For the formulation of the following results we need to introduce a new set of parameters. 



x^ Pnix)W (x) dx , k = n,n + l,... 



(2.20) 



where Pn are the monic MOPRL. In the spirit of [48], we can express them in terms of the 
coefficients a„j- of the polynomials Pn (see (jl.ip ) in a form suitable for numerical implementa- 
tion: 

Lemma 2.10. Let CI be the block lower triangular matrix built up from the coefficients of the 

MOPRL {Pn)n, 



( I 



so that 



fPo{x)\ 
Pi{x) 

\Pnix)J 



( I \ 

xl 
\x'^l) 



I = In. 



Then b^^k defined in (|2.20p can be obtained from the last N x N block row of ft ^ as follows: 

K^Knln-k = k = 0,...,n. (2.21) 
Proof. Observe that 

(^0,n) • • • ) ^n-l,n? ^n/n) = / x"'W (^Pq , Pi , . . . , P*) dx 

J a 

= I W (I, X/,..., X"I) dx = (/X„, //„+!,..., //2n)r2*, 

J a 



Properties of Matrix Orthogonal Polynomials via their Riemann-Hilbert Characterization 9 



where fj-n = J^^ x"'W{x)dx are the moments of the weight matrix W. Analogously, by (|1.2p . 



(diag (7fc)Lo)"^ = diag {%Yk=o = / 

J a 



fPo\ 

Pi 



w{pip,*,...,p:)dx 



/ I \ 

xl 



\x"lj 



w(i,xi,...,x''i)n* = n 



n*. 



Hence, 



(^0,n70, • • • , K-l,n'yn-l, K.,nln) 

Ml 



-1 



n-i = (o,o,...,j)n-\ 



\A*n Mn+l • • • ^J'2n / 



which yields ifMT]) . 



For what follows it is useful to single out the explicit expressions for bn~3,n, ■ ■ ■ , bn,n, that 
are obtained from Lemma 12.101 by direct computations: 

Corollary 2.11. Coefficients bn-k,n> for k = 0, 1,2,3, are given by 

'yn-lOn-l,n — ~0-n,n-l^ 

7n-2&n-2,n = 0.*n-l,n-2^*n,n-l ~ ^*n,n-2i 

In-Z^n-Z^n = — a.,i_2,n-3^n-l,r!.-2^n,n-l + ^n-l,7i-3*^n,n-l + *^n-2,n-3^n,7i-2 ~ *^n,n-3- 

Now we return to Theorem 12. 2t as its immediate consequence we can relate the coefficients 
in the asymptotic expansion of and of {Y^^ ^ ^ with the coefficients o,nj and ^nj* 



Corollary 2.12. The coefficients in (12. Sh are given by 

1 . 



V." 



ar 



n ."n,n.+i— 1 



i > 0. 



(2.22) 



-27ri7„_ia„„i_„,_ 
Analogously, the coefficients Y-^ in (|2.5p are given by 



( 



i > 0. 



(2.23) 



\^ 27ria„_]^ ; ci„ 



Remark 2.13. One of the consequences of (|2.14p is that a„j- and bnj appearing in ()2.22p . ()2.23p 
belong to M^^^. 

The explicit expressions for coefficients of the asymptotic expansion above, combined with 
the obvious fact that Y'^{Y")~^ = I2N, yield in a straightforward way the following identities: 
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Proposition 2.14. The coefficients of the monic MOPRL Pn and the coefficients of the Cauchy 
transform C{PnW) {i.e. the sequences bn^k) satisfy the following relations 

m m 

an,n-m+j^n-l.n+.7-l = bn^n+m~j~lO'n-l,n-jJ m > 1, (2.24) 

j=0 j=0 

and 

m m 

^ ^ 0'n,n—m+jbn^n+j-l ~ ^ ' ' _, ^n,n+m—j—lO-n,n—j^ m > 1. (2.25) 

i=0 j=0 

In particular, for m = 1 we recover the second identity of Corollary 12.111 
Proof. From the asymptotic behavior of 1^" and (Y^)~^ and using Y^{Y'^)~^ = I we get 

m 

^Y^_^Yp = 0, m>l, 

3=0 

from which block entries (1, 1) and (1,2) give ()2.24p and ()2.25p . respectively. ■ 

We finish this section with the matrix analogue of the Christoffel-Darboux (CD) formula for 
MOPRL, written in terms of the solution of the characterizing RHP. For an orthonormal family 
{Pn)n we consider the kernel 

n-l 

Kn{x,y) = Y,Pj{y)Pjix), x,yeR. (2.26) 

j=0 

Then the CD formula (see 12011 reads as 



, P* Jy)AnPnix) - P*(y)A*Pn^i(x) 

Kn{x,y) = ""^^^ ^ ^ n n ) ^ ^.27 

x-y 

Proposition 2.15. The matrix CD kernel (|2.26p satisfies, for x,y ^R, 

Kn{x,y) = -—^ -iYiM'YiAx) - Y^M^2i{x)) 

2Tn[x — y) 



2-Ki{x 

The proof follows taking into account the Christoffel-Darboux formula (]2.27p . the expression 
for An in (j2.17p and the expression of the inverse (1^")~^ in (j2.16p . The last proposition for 
general weight matrices of rectangular size can be found in [18j . 

2.2 Difference and differential relations 

Both the MOPRL and the solution of the RHP given in Theorem 12.21 depend on two variables, 
the discrete n and the continuous z. We can derive further properties by analyzing the variation 
of this solution with respect to either variable, which yields linear relations of the form 

Y'^-^\z) = En{z)Y^{z), ^^"(^) = Fn{z)Y^{z). (2.28) 

The general methodology to get these equations may be traced back to the original work of 
Gelfand, Levitan and other authors. The idea is that if the jump matrix for a RHP is independent 
of a variable, then a variation with respect to that variable leads to an identity. The fact 
that the jump matrix in (j2.2p is independent of n allows one to obtain immediately the first 
identity in (I2.28p . which in turn is connected to the well-known three-term recurrence relation 
for MOPRL. 
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Theorem 2.16. There is a unique sequence of matrix coefficients cxn such that the solution of 
the RHP satisfies the following first-order difference equation 

= ( ^[~^;^\2^^l_ I Y'^iz). (2.29) 
-27rf7„ i 

Consequently, the monic MOPRL P„ satisfy the following three-term recurrence relation: 

xPnix) = Pn+i{x) + a„P„(x) + /3„P„_i(x), n > 0. (2.30) 

The coefficients q;„ and as well as -jn {defined in (j2.7p ) can be expressed either in terms of 
the elements of the solution of the RHP or in terms of the coefficients (jl.ip ) as follows: 

Proof. The matrix-valued function R = >^"+i(>^")-i^ analytic in C \ [a, 6], satisfies i?+(x) = 
R-{x) for all X € {a ,b). By the Morera's theorem, R is analytic in C \ {a, 6}. The behavior at 
the exceptional points gives that the singularities at a and b are removable, and R is an entire 
function. From (12.31) and Liouville's theorem we conclude that 



-n+lC^N _ ' V^l III V'l III 

" " ' (n-^'); 




y"(z). (2.32) 



This proves (j2.29p along with the expression for «„, in terms of the solution of the RHP. The three 
term recurrence relation ()2.30p for the MOPRL is obtained by considering the (1,1) block entry 
of (j2.32p . Finally, the rest of the expressions for the coefficients is found using Corollarv 12.121 
(formula ^^). ■ 

We turn now to the dependence of Y"^ on the continuous variable z; notice however that the 
jump ()2.2p along the real line does depend on this variable, so the application of the paradigm 
explained at the beginning of this subsection is not straightforward. At this stage we can only 
aspire to find a differential relation like in (j2.28p . but with a matrix depending upon the 
entries of 1^" (or equivalently, upon the MOPRL themselves). 

Theorem 2.17. The solution of the RHP for Y'"'{z) satisfies the following first-order matrix 
differential equation 

r"(z) = I i.~2v^7^' I ^"(^). (2-33) 



dz 
where 



27ria„_i(z)T„-i; 



t=a . 



and 



rb 



Pn{t)W'{t)P*_,{t) ,^ P^{t)W{t)P*_,it) 



^n{z) = -l '''' "-^^^ dt 



t=b^ 



t=a . 



(2.34) 



7n-i, (2.35) 



provided that all the functions in the right-hand sides of (|2.34p and (j2.35p exist for z € C \ [a, 6]. 
In particular, 

7n2i;(z) = %n{z)ln- (2.36) 
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Remark 2.18. Formulas (j2.34p and (|2.35p were introduced in the scalar case in [7] and in 
the matrix case (for a normalized family) in [23]. In general, they bear a formal character. 
In particular, the boundary terms in (j2.34p and (j2.35p may be undefined, depending on the 
behavior of W at the endpoints (see (j2.ip ). As it will be shown in the next section, these 
formulas simplify considerably after considering a special transformation of the RHP for Y^. 
Finally observe that the coefficients (j2.34p and (j2.35p are unique since Y"^ is unique. 

Proof. Consider the matrix-valued function RP'{z) = [^^"(-2)] From the explicit 

expressions (j2.15p and (j2.16p we have that the block entries of are 

(fl"(z))n = -27ri(P^(z)C(\yP;_i)(^) -C(P„Vr)'(z)P:_i(z))7„_i, 
(fl"(z))i2 = -PL{z)C{WP:){z) +C{PnW)'{z)P:{z), 

(fl"(z))2i = (2vri)V-i(^n-i(^)C(WP,:-i)(^) -C(P„_iVr)'(z)P;_i(z))7n-i, 
(l?"(z))22 = -2vri7„_i( - P:,_^{z)C{WP:){z) + C{Pn-iW)' {z)P:{z)) . 

In particular, -(i?"(z))ii = (l?"(z))^2 and -27ri7„(P"(z))i2 = ^^{R^+^{z))2il-y 

We will now derive formulas ()2.34p and (j2.35p . For that purpose we need the following 
technical observations: 

Lemma 2.19. 

(i) For every P € 

P{z)C{WP:){z) = C{PWP:){z), C{PnW){z)P{z) = C{PnWP){z). 
(a) For any differentiable and integrahle matrix function F bounded at z = a, b, 



dz 27ri t — z 



z G C \ [a,6]. 



Identities in (i) follow by adding and subtracting C{PWP*){z) and using the orthogonality 
of the MOPRL, and (ii) is obtained by integration by parts and using that ^(l/(i — z)) = 

We return to the proof of Theorem 12.171 showing in detail how to obtain 5S„(z). For 2l„(z) 
the computations are similar and will be omitted for the sake of brevity. 

Using (i) from the previous lemma and the product rule of differentiation we get 

(i?-(z))n = -2vri(C(P^l^P;_i)(z) -C(P„W-P„*_i)'(z) +C(P„l^(P„*_i)')(z))7n-i. 
Applying now (ii) to the second term and canceling we get 

1 Pnit)W{t)P*_,{t) 




(P"(z))n = -2vri \^-C{PnW' P:_^,){z) + ^ _ ^ 

which yields ([2:35]) . ■ 

Let us discuss now two main consequences of the existence of the first order matrix-valued 
differential equation ()2.33p . 

One one hand, equations (j2.28p . known as the Lax pair for the RHP, are clearly overdeter- 
mined, so compatibility conditions (via cross-differentiation of both equations) yield 

E'^{z) + En{z)Fn{z) = Fn+i{z)En{z), (2.37) 
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also known as string equations. They allow us, for instance, to recover the sequence {F^} 
from the sequence {En} and the initial value -Fi. These ideas relate also the RHP to nonlinear 
problems such as the nonlinear Schrodinger equation (NLS), see [T7j, and the Toda flow, see j27j . 
In our situation, the compatibility conditions (|2.37p considered entry-wise imply the following: 

Proposition 2.20. The recurrence relations 

I + *B„+i(z)(zI - a„) - (zl - a„)5S„(z) = 2i:+i(z)/3„+i - PnK~i{z) (2.38) 

and 

^n+l{z) + 7n 'K(^)7n = {zl " CXn)K{z) (2.39) 

hold for every n > 0, where q;„ and Pn oif^ the coefficients of the three-term recurrence rela- 
tion (fOOl) . 

Proof. Block entries (1, 1) and (1,2) of (|237D give (|238D and (f239D . respectively. Block entry 
(2, 1) of (|2.37p is equivalent to (|2.39p using relation (j2.36p and q;„ = 7~^q;*7„ (which is a con- 
sequence of substituting ()2.22p for z = 2 in ()2.29p as z ^ oo, along with (j2.3ip ). while block 
entry (2, 2) is equivalent to (|2.38p as well. ■ 

Another consequence of Theorem 12.171 is the existence of the ladder operators for MOPRL: 

Corollary 2.21. The monic MOPRL {Pn)n satisfy the following difference- differential relations 
{lowering and raising operators, respectively): 

Pl,{z) = -^n{z)Pn{z) + Kiz)(3nPn~l{z) (2.40) 

and 

Pni^) = [K{^){ZI - OCn) - ^n{z)]Pn{z) - K{z)Pn+l{z). (2.41) 

Proof. Block entry (1, 1) of ()2.33p gives the lowering operator ()2.40p . using the expression for /3„ 
in Theorem 12. 16t while block entry (2, 1) of (|2.33p gives the raising operator (|2.4ip using (|2.36p 
and (f239D . ■ 

Proposition 12.201 and Corollary 12.211 were already known for the normalized MOPRL in [J^ ; 
the RH approach gives an alternative proof of these results. The ladder operators are the basic 
differential relations for MOPRL. It is well-known that they can be combined to build a second- 
order differential equation satisfied by the polynomials. This fact was mentioned in [25], but no 
explicit expression of the differential equation was given; here we present it for completeness: 

Corollary 2.22. The polynomials {Pn)n satisfy the following second-order differential equation 

P'^{z) + Tln{z)P^{z) + ^n{z)Pn{z) = 0, (2.42) 

where 

mn{z) = -(2t;(z))'2t-*(z) + 5S„(z) - K{z){zl - a„) + 2i;(z)*B„+i(z)2t;:*(z) 

and 

9I„(z) = 5Jt„(z)5S„(z) - ^l{z) + Kiz) + %{z)(3nK-i{^), 

provided that the inverse o/2l*(z) exists for z £ C\[a,b]. Here and below, the superscript —* 
denotes the conjugate transpose of the inverse. 
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Proof. Differentiate the lowering operator ()2.40p and substitute the raising operator ()2.4ip 



Remark 2.23. Note that we can easily obtain another second-order differential operator sat- 
isfied by the orthogonal polynomials reversing the ladder operators. This new differential equa- 
tion needs not in principle be the same as (j2.42p . Nevertheless, it is straightforward to see that 
both equations are equivalent using the compatibility conditions ()2.38p and ()2.39p . 

Although Theorem 12.171 gives explicit expressions of 2t„(2;) and ^niz), these coefficients are 
usually difficult to calculate. They are not even defined on the interval [a,b]. In the next 
section we will introduce some additional assumptions that simplify the differential equation 
considerably. 

3 Transformation of the RHP when supp(VF) = R 

As it was pointed out in the previous section, the independence of the jump in (j2.2p with respect 
to the discrete variable n allows one to obtain the three-term recurrence relation in a straight- 
forward way. Under additional assumptions on the jumps we can perform a transformation of 
the RHP in such a way that a new jump is independent of the continuous variable z. This will 
have consequences on the resulting differential relations. 

In the matrix case there is some extra freedom absent in the scalar situation, that gives us 
a whole new family of differential relations, and consequently, a whole class of ladder operators, 
some of them of the 0-th order, something that is not possible in the scalar case. The combination 
of the ladder operators will give rise to a class of second-order differential equations, some of 
them of order less than 2. 

For simplicity, we will focus here on the case when [a, b] = M, assuming along this section 
that W is smooth and positive definite on the whole real axis M. Some ideas about how to handle 
the case of finite endpoints of the support of W are briefly explained in the final Section [5l 

3.1 The transformation 

As we have seen in Section [2.21 the recurrence relation (first identity in (|2.28p ) is a general fact 
intrinsic to the orthogonality of the polynomials. Simple differential relations are much more 
demanding, and we must impose at this stage further conditions on the orthogonality weight W. 

Our immediate goal is to obtain an invertible transformation such that X" has 

a constant jump across M. We will consider 



where V \s & matrix-valued function, analytic in C \ M and continuous up to M, and invertible 
for all z € C. Then the jump matrix for X" on R is 



Observe that this kind of transformations does not affect the first difference equation in the 
Lax pair (j2.28p . but constant jumps will allow us to use the strategy of variation of the problem 
along z. Going down the path of simplification, we can try a block-diagonal matrix 



evaluated at n — 1. 




(3.1) 
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where T is an invertible N x N matrix function. Then (j3.ip boils down to 
IT^WT 

These preliminary considerations motivate to consider a factorization of the weight in the 
form 

W{x) =T{x)T*{x), xeM, (3.2) 

where T is a smooth matrix-valued function on M. Under our assumptions the existence of such 
a T is guaranteed, but not its uniqueness: performing its RQ decomposition it can be written as 

T{x) =f{x)S{x), xGR, (3.3) 

where T{x) is an upper triangular matrix and S{x) is an arbitrary smooth and unitary matrix 
(for each x S M). This representation can go further taking into account that any unitary matrix 
can be written as S{x) = e^^^\ where Q(x) is an skew-Hermitian matrix function. Additionally, 
since any skew-Hermitian matrix is normal, it has a factorization Q{x) = iU{x)D{x)U* (x), 
where U is again unitary and D a diagonal matrix with real entries. Therefore T from (|3.3p 
can be written as 

T(x) = f(x)C/(x)e^'°(^)j7*(2;). 

We narrow the choice of T (and W) by imposing the additional constraint that there exists 
a matrix polynomial G such that in a neighborhood of the origin, 

T'{z) = G{z)T{z). (3.4) 

As a consequence, (see the discussion in [39l Chapter 9] or [49l Chapter 1]), T has an analytic 
continuation to the whole plane as an entire and invertible matrix-valued function, and 

W{z) = T{z)T*{z) 

provides an analytic extension of W to the whole plane. 

Following our previous discussion, we define the entire matrix-valued function 



V{z) 



T{z)\ 0. 

:T"*f 



Observe that V is invertible for all z € C. By differentiating the identity T = I and 
using dM]) we conclude that (T-^)' = -T'^G, so that 

, ( G{z)\ \ , , , ^ 

^ = ( o'i-Gni)-] 

If Y^{z) is the solution of the original RHP, let us define 

X"(z) = Y"(z)y(z). (3.6) 

Then straightforward computations show that X''^{z) is analytic in C \ M, satisfies the jump 
condition 

X-(x) = X!!(x)f-^.i^^-V xGM, 



:/ 



N 
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and for m G N (see ([O]) ). 



.m+l \ 



i=l 



6 







Viz), 



Z ^ CO. 



(3.7) 



Since X"(z) is invertible in C\IR, we can consider the matrix function Fn{z)= [^X"(z)]X"(z)~^. 
Again, Fn{z) is analytic in C \]R and on the real line, (i^„)+(a;) = (i^„)„(a;), which implies that 
it is an entire matrix function. From (j2.3p and (j3.7p . for z — )• oo. 



dz 



[Y-{z)r' = o{^^ 



and combining it with (|2.3p . (j2.5p . (|3.5p and ()3.7p . we get for z — > cx). 



i^n(^) 



X"(z)-^ 



+E 1^ + ^(1/^ 



m+l^ 



G{z) 








-G*(z) 



^27V + E ^ + 0{l/z-^+^) + 0{l/z). 



(3. 



By Liouville's theorem, the right hand side in (j3.8p will coincide with its polynomial part in the 
expansion at infinity, and its degree is no greater than the degree of G. To be more precise, if 
we assume that the degree of G is m G N and denote 



Mi 











G{z) = ^MjZ^, My. 

j=0 

then after dropping the negative powers of z in 

i^n(.)=(^l2^ + 5:^ + 0(l/.-+^)j 



M* j ' 



,3=0 



k=l 



we obtain that 



m / m j—k 

Fn{z) = j2[Y.j2'^-M,Yr.,., 

k=0 \j=k i=0 



with l^o" = l^o" = I27V. 



For instance, if m = 0, 

Fn{z) = Mo, 
for m = 1, 

F^{z) = Miz + Mo + l"i"Mi + MiYi", 
and for m = 2, 

Fn{z) = M2Z^ + (Ml + Yi^Ma + M2Y{')z + Mq + Yi^Mi + MiY^" 
+ 1^2" ^2 + Y{'M2Y{' + M2l^2''- 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



Finally, using the explicit expressions (j2.22p and (j2.23p in (j3.9p . we arrive at the following 
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Theorem 3.1. Under assumptions (j3.2p and (j3.4p . with G{z) = ^ MjZ^ G Pm, i/ie matrix 

j=0 

function JC" defined in (I3.6p satisfies the following first-order differential equation with polyno- 
mial coefficients: 



d_ 

dz 



X^{z) = F^{z-G)X^{z), 



where 



Fn{z-G) 



-Bn{z-G) 

2TTi An^iiz; G)jn~l 
An and Bn are matrix polynomials, 



-— 7„M„(z;G) 



(3.13) 



(3.14) 



B:{z;G) 



' m—l 



An{z;G) = -7„ I bn,n+m-j-l^lniz) + ^j,n{^)b*^,n+m,-j-l ) > 

,J=0 



+m-j-2A*. 



(3.15) 
(3.16) 



and ifte coefficients A.j^niz) are given by 



^J>(z) = Y Pn,k{^)Mjn-j+k, Pn,k{z) = Z^ I + an,n-lZ^ ^ "I h a„,„_ 



fc=0 



Moreover, Fn{-;G) is linear in G: 

Fn{z; Gi + G2) = Fn{z- Gi) + Fn{z; G2) 



(3.17) 



Remark 3.2. The coefficients bn^k were introduced in (|2.20p and discussed in Lemma [2. 101 It is 
worth observing the similarity of this result with Theorem I2.17| in the present situation we can 
compute the differential equation for X"" directly in terms of the coefficients of the MOPRL Pn 
without considering integrals or studying the behavior at the endpoints. Finally, we have once 
again that 



7„A* (2;; G) = An{z; G)7„. 



(3.18) 



Again, for lowest degrees m in Theorem 13.11 we can use (jS.lOp . (|3.1ip . (j3.12p . Lemma 12.101 
and Corollary 12.111 in order to write the coefficient matrix p.l4p of the differential equation 
explicitly. For instance, for m = 0, 



Fn{z;G) 
for m = 1, 



Mo 

6 







Fn{z;G) 



1 
27ri 



(7-^Mi*+Mi7-^) 



-G*{z) + al^^_,M;-M^a 



n,n— 1 



G{z) + an,n-lMi - Mian,n-l 

'-2nii^n-iMiTM[-^n-^''' 
and for m = 2, with the notation ()3.14p . 

Bn{z; G) = -G{z) - (a„ ,„_iM2 - M2a,„,„_i)z + Mia„ - a„ ,„_iMi - a„ ,„_2M2 



, (3.19) 
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An{z; G) = -M^z - Ml + <+i^„M2* - MK^^.^ 

- 7„(M22 + Ml - M2an+l,n + an,n-lM2)j-\ 

Observe from Theorem 13.11 that Bn{z;G) is a matrix polynomial of degree at most m 
and An{z\ G) is a matrix polynomial of degree at most m — 1. 

Finally, it is important to emphasize that in many situations we can exploit the implications 
of the freedom in the factorization (j3.3|) on the differential equation (|3.13|) . This freedom, as we 
will see in Proposition 13.91 appears only in the matrix setting. 



Proposition 3.3. If under assumptions (j3.2p and p.4p . with G a polynomial, there exists 
a non-trivial matrix-valued function S, non-singular on C, smooth and unitary on M, such that 

H{z) =T{z)S'{z)S*{z)T~\z) (3.20) 

is also a polynomial, then T = TS satisfies 

W{x) = T{x)T*{x), X G M, T'{z) = G{z)T{z), z G C, 

with G{z) = G{z) + H{z). Moreover, the matrix X" defined in (j3.6p . satisfies along with 
(j3.13p - (j3.16p the following relation: 



^X^{z) = {Fn{z; G) + Fn{z; H)) X^iz) - 



with x{z) = S'{z)S*{z). 














(3.21) 



Remark 3.4. Equations (j3.13p and (|3.2ip are not necessarily trivially related, and in principle 
we could combine them in order to establish new relations for X"^, and thus, for Y^. 

Proof. The key observation is the formula (|3.17p . so that the coefficients An and Bn in (|3.13p 
depend linearly on G: 

An{-;G + H) = An{-;G) + An{-,H), Bn{-,G + H) = Bn{-,G) + Bn{-,H). (3.22) 



3.2 Differential properties 

Taking into account the similarities between Theorems 12.171 and 13.1} we get immediately the 
analogue of Proposition 12.201 (the compatibility conditions): 

Proposition 3.5. Under assumptions (j3.2p and (|3.4p . with G{z) a polynomial, the coefficients 
of the differential equation (j3.13p satisfy the following recurrence relations: for every n > 0, 

I + Bn+iiz; G){zl - an) - {zl - cXn)Bn{z- G) = ^:+i(z; G)f3n+i - PnAl_^{z- G) (3.23) 

and 

Bn+i{z; G) + in^Bliz; G)7n = {zl - cxn)A*n{z; G), (3.24) 
where an and (3n o,Te the coefficients of the three-term recurrence relation (j2.30p . 

As before, the ladder operators (the most basic differential properties for MOPRL) can be 
easily obtained by analyzing the first block column of X^ in the differential equation (|3.13p : 
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Proposition 3.6. Under assumptions ()3.2p and (j3.4p . with G{z) a matrix polynomial, the 
monic MOPRL {Pn)n satisfy the following difference- differential relations {lowering and raising 
operators, respectively): 

p:,{z) + Pn{z)G{z) = -Bn{z; G)Pn{z) + Al{z; G)f3nPn-l{z) (3.25) 

and 

P^{z) + Pn{z)G{z) = {Al{z-G){zl - cxn) - Bn{z;G))Pn{z) - Al{z-G)Pn+i{z), (3.26) 
where An and Bn are given by (|3.15p and (|3.16p , respectively. 

Proof. Block entry (1, 1) of (|3.13p gives the lowering operator (j3.25p . using (|2.3ip and (|3.18p . 
while block entry (2, 1) of ()3.13p gives the raising operator ()3.26p using (|3.24p and ()3.18p . ■ 

Remark 3.7. Comparing with the results of Corollarv l2.2H notice that these ladder operators 
contain a term with the MOPRL multiplied on the left. 

In the situation described in Proposition l3.3| we can in principle exploit the non-uniqueness of 
the relations above in order to derive further relations for the MOPRL. The following result shows 
the possibility of the existence of 2-terms recurrence relations for the family Pn, a phenomenon 
that has not been reported before in the theory of MOPRL. 

Corollary 3.8. Under conditions of Proposition [3^ the family of monic MOPRL {Pn)n satis- 
fies, along with (I3.25|) - f3.26p . the following relations: 

Pn{z)H{z) = -Bn{z; H)Pn{z) + ^^(z; H)(5nPn-l{z) (3.27) 

and 

Pn{z)H{z) = {Al{z-H){zl - CXn) - Bn{z-,H))Pn{z) - Al{z-H)Pn+l{z), (3.28) 

where H is defined in ()3.20p . 

In particular, provided An{z; G) is invertible for all z £ C, 

Pn{z)H{z)+Bn{z;H)Pn{z)-Al{z-H)An*{z-G) 

X (P^(z) + Pn{z)G{z) + Bn{z; G)Pn{z)) = 0. (3.29) 

Proof. To get (|3.27p . we subtract (j3.13p and (|3.2ip and then evaluate the (1,1) block entry. 
(j3.28p is obtained from (I3.27P and the three term recurrence relation. Finally, (I3.29P is obtained 
by replacing ([3:25]) in ([3:271) . ■ 

Equations (|3.27p and ([3.28p are known as the 0-th order ladder operators, while ([3.29P is 
a first order differential relation for the MOPRL. In some situations they yield trivial identities; 
this is always true in the scalar case, as the following proposition shows: 

Proposition 3.9. Assume that under the conditions of Proposition^^ x{z) = S' {z)S*{z) = 
ip{z)I, where p is a scalar polynomial of degree m. Then 

An{z;H) = 0, and Bn{z; H) = —ip{z)I . 

Proof. From ([3:20]) . H = TxT"^, and by hnearity oi An and Bn in H (see ([3:22]) 1. it is enough 
to prove the formulas above for monomials z^'I. In this case, formulas (|3.15p and (|3.16p simplify 
considerably. Using Proposition 12.141 gives An{z; H) = and Bn{z;H) = —iz^I. ■ 
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In the scalar case {N = 1) the only smooth and unitary function s{x) on M has the form 
s{z) = e^P^^\ with p real-valued, and 

s'{z)s{z) = ip{z), 

so that the assumption that p is a polynomial brings us to the situation described in Propo- 
sition 13.91 This explains why in the scalar setting we never get nontrivial 0-th order ladder 
operators. In the scalar case one cannot have first order differential relations for the MOPRL, 
such as those given in (j3.29p . The general matrix case is much richer and complex, as will be 
illustrated with some examples in the next section. 

We finish this section with a class of second-order differential equations satisfied by the 
MOPRL (-Pn)n- As a consequence of the Proposition 13.61 we have the following 

Proposition 3.10. Under assumptions ()3.2p and (|3.4p . with G{z) a polynomial, the MOPRL 
{Pn)n satisfy the following second-order differential equation 

P'^ + 2P^G + Pn [G' + G^) + MnPL + MnPn + MnPnG = 0, (3.30) 

where An{z) = An{z; G), Bn{z) = Bniz; G), 

Mn{z) = Mn{z;G) 

= -{Al{z))'A~*{z) + Bn{z) - AUz){zI - CXn) + A*^{z)Bn+i{z)A~* (z) , 

and 

Mn{z) = Mn{z; G) = Mn{z)Bn{z) - B^z) + B'^{z) + Al{z) (3nAl^,{z) , 

provided that the inverse of An{z;G) exists for z € C. Here, again, —* denotes the conjugate 
transpose of the inverse. 

Proof. Differentiate the lowering operator ()3.25p and substitute the raising operator ()3.26p 
evaluated at n — 1. ■ 

Remark 3.11. Note that we can easily obtain another second-order differential operator sat- 
isfied by the orthogonal polynomials reversing the ladder operators. This new differential equa- 
tion needs not be in principle the same as (j3.30p . Nevertheless, it is straightforward to see that 
both equations are equivalent using the compatibility conditions ()3.23p and ()3.24p . 

The equation (|3.30p does not have the form of the right hand side differential operator 
considered for instance in [23], due to the terms MnPn, J^nPn and MnPnG. In some cases, 
under additional assumptions on the weight, ()3.30p can be reduced further, as we will see in the 
following section. 

4 Illustrative examples 

In this section we study a number of examples of weights W which are smooth and non-vanishing 
on the whole real line; this assumption simplifies the Riemann-Hilbert formulation because in 
this case one does not consider the local conditions (Y4) (see Section 12. ip . Additionally, without 
loss of generality, we take VF(0) = /. 

For convenience, we consider the weights of the form 



Vr(x) = e-23(^)l7(x)f7*(x), xGM, 



(4.1) 
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where g is a scalar real- valued function, so that with the notation (j3.2p and (|3.4p . we may take 
T(x) = e-^(^')j7(a;) and G{x) = -q'{x)I + U'{x)U-^{x). (4.2) 

We are interested in the case when G is a matrix polynomial. Hence, keeping up with previous 
hypotheses, we will assume that g is a (scalar) polynomial of even degree with real coefficients 
and a positive leading coefficient, and that U'{x)U~^{x) is a matrix polynomial. 

Since our main goal here is to generate a set of examples (some new, some already well known), 
we restrict the degree of q to either 2 (the Hermite case) or 4 (the Freud case) , and the degree of 
U'{x)U~^{x) to at most 1. We start by considering the case when U' {x)U~^ (x) is a monomial, 
i.e. either U'{x)U~^{x) = A or U'{x)U~^{x) = 2Bx for constant matrices A, S G C^^^. 
Taking into account the linearity ()3.22p of the coefficients An and Bn in (jS.lSp . (I3.16p . we 
can obtain the differential equation (j3.13p for the general case of U'{x)U~^{x) = A + 2Bx. 
However, finding the corresponding orthogonality weight is more involved: when A and B do 
not commute, solving U'{x)U~^{x) = A + 2Bx is not straightforward. In Section [4.1.3l we will 
discuss some examples, which yield explicit expressions for U{x), related to this case. 

Recently an example has been found in [3] where a weight matrix supported in the real line 
is explicitly given (but not of the type (14. ip ). when G{x) is a matrix polynomial of degree N 
with in general non-commuting coefficients. 



4.1 The Hermite case 

For q{x) = x^/2, let us consider two cases, first U' {x)U~^ (x) = A and then U'{x)U~^{x) = 
2Bx. We end up this Section by discussing briefly the case of U'{x)U~^{x) = A + 2Bx. 

4.1.1 U'{x)U-^{x) = A 

The differential equation U'{x)U~^{x) = A, so that U{x) = e"^^, T{x) = e~^^/^e"^^, and the 
weight matrix (j4.ip is given by 



W(x) = e-^'e^^e^*^ A G C^^^, x G M. (4.3) 
The matrix G has the form 

G{x) = -xl + A, 
and according to ()3.19p . 

Anix; G) = 21, Bn{x; G) = -G{x). 
The compatibility conditions of Proposition 13.51 yield 

2(/3„+i -(3n) = I + Acxn - anA, (4.4) 

and 

a„ = ^(A + 7^U*7„). (4.5) 
The lowering and raising operators from Proposition 13.61 are reduced now to 

P^(x) + Pn{x)A - APn{x) = 2/3„P„_i(x), (4.6) 

and 



-P^(x) + 2xPn{x) + APn{x) - Pn{x)A - 2«„P„(x) = 2P„+i(x). 
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Summing up the telescopic relation in (j4.4p (or comparing the 0{x'^ ) term in ()4.6p ) gives 

f3n = ^{nl + an,n-iA - Aan,n~l)- (4.7) 

With the notation of Proposition I3.1UI 

Mn{x) = 2 (a„ - A) , J^n{x) = -2 (a„ - A) G{x) - G\x) +1 + 4/3„, 

so that the differential equation (|3.30p for the monic polynomials P„, orthogonal with respect 
to the weight (j4.3p . boils down to 



P^{x) + 2P^(x)(A - xl) + Pn{x){A^ - 2xA) 

= {-2xA + A2 - Apn)Pn{x) + 2(A - an)(Pn(x) + P„(x)A - APnix)). (4.8) 

These formulas hold for any constant matrix A, and we cannot expect important sim- 
plifications without narrowing the class of the weights further. This can be done assuming 
in addition that the hypotheses of Proposition 13.31 hold. This, as it was shown in \23\ 124] . 
imposes additional constraints on the weight W. Since the construction is described in detail 
in \23\ I24j . the exposition in this part will be rather sketchy. 

The matrix H from (13.20p . restricted to R, can be written as 



//(x) = T(x)xT-Hx) = e^^Xe^ 



Ax 



= X + ad^(x)x + adi(x)y + • • • = E ""^Aix)^- (4-9) 

Here x(x) = S'{x)S*{x) is skew-Hermitian on M, ad^ is the commutator given by adA{x) = 
Ax — xA, and we define recursively 

^d^Aix)=X, ^d'^tHx) = adAiad^ix)) for n > 1. 

The simplest situation obtains when the right hand side in (|4.9p is constant; as it was shown 
in Lemma 2.4 of [24], this assumption yields x = ^^-^ for certain a € M (which is the case 
discussed in Proposition 13. 9p . when there are no new ladder operators. Consequently, for the 
first non-trivial situation we must assume that ()4.9p is a matrix polynomial of degree at least 
one. We consider two situations that yield degree exactly one (see [23\ [23] for motivations and 
further details). 

Let us define a nilpotent matrix of the form 

N-l 

L=Y^ VkEk,k+i, i^fc € C \ {0}, (4.10) 
fe=i 

where Eij is a matrix with 1 at entry (^,j) and elsewhere, and a diagonal matrix 

N 

J = J2{N-k)Ek,k. (4.11) 

k=l 

For the first non-trivial example we assume that A = L and x = '^Ji so that 

adA(x) = -A (4.12) 
and ad\{x) = 0. Thus, H{x) = e^^xe""*""' = i{J - Ax), and by ([3T9]l . 
Anix; H) = i{-A* + 7nA7-^) = 2i« - A*), 
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13n{x; H) = i{Ax - J + a„_„_iA - Aa„,„„i) = i{Ax - J + 2(3n - nl). 
From the compatibility conditions of Proposition 13.51 we get 

Ja„ - a„ J + a„ = A + 1 ( A^a. - a„ A^) and 

J - 7^^ J7„ = Aq!„ + q;„A - ^ (4.13) 

where we have used relations (|4.4|) . (j4.5p . (|4.7|) and 

(A - 2q!„)/3„ = /3„(A - 2q;„_i), 

which can be easily proved using (|2.31|) and (|4.5|) . 

Thus, from Corollary 13.81 (j4.4p and (j4.13p . the lowering and raising operators (of the 0-th 
order) are 

Pn{x)J - JPn{x) - x{Pn{x)A - AP„(x)) + 2/3„P„(x) - nP„(x) 
= 2(A-a„)/3„P„_i(x), 

and 

Pn{x){J - xA) - 7-i(J - xA*)7„P„(x) + 2/3„+iP„(2;) - {n + l)P„(x) 

= 2{cXn - A)Pn+l{x), 

respectively, which yields by (|3.29p the first order relation 

(A - cxn)PrXx) + (A - cxn + x/)(P„(x)A - AP„(x)) - 2/3„P„(x) 

= P„(x)J- JP„(x)-nP„(x). (4.14) 

These identities hold in addition to the second order equation obtained in (|4.8p . valid as we 
recall, for any A. As far as we are aware of, these relations are new. We use them to simplify 
the differential equation remarkably: multiplying ()4.14p by 2 and plugging it into (j4.8p gives us 
for A = L, 

P^'(x) + 2P^(x)(A - xl) + Pn{x){A^ - 2J) = (-2n/ + A^ - 2j)P„(x), 

which is a linear second-order differential equation with coefficients in the right hand side 
independent on n (a.k.a. Sturm-Liouville equation with polynomial coefficients), studied by 
Duran and Griinbaum in [23] . 
In the second case we take 

N-l 

A = L{I + L)-^=Y,{-iy~^L\ x = iJ: 
i=i 

where L is given in ()4.10p and J in ()4.1ip . Consequently, 

adA(x) = -A + A\ adi(x) = 0, (4.15) 

and 

H{x) = e^^xe"^^ = z(j - (A - A^)x). 
In this case, 

An{x; H) = i{-A* + {A*f + 7n(A - A^)-^-^) 
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= 2i« - A* - « - A*X - - A*)) 

and 

Brr{x;H)=i{{A-A^)x-J + ar,,n-i{A-A^) - (A- A2)a„,„_i) 
= i((A - A2)x - J + 2/3„ - nl - 2(AA^ + /9nA) + 2nA). 

As to new results, we get the lowering operator, 

Pnix){J - x{A - A2)) - (J - x{A - A2))P„(x) 
= i-2(3n + n{I - 2A) + 2(A/3„ + /3„A))P„(x) 

- 2(q„ - A-{cXn- A)OLn - OLn{oLn " A))/3nP„_i (x) , 

the raising operator 

P„(x)(j - x{A - A^)) - (J - x(A - A2))P„(x) 
= (-2/3„ + n{I - 2A) + 2(A/3„ + /3„A))P„(x) 

- 2(q„ - a - (q!„ - A)q!„ - Q!„(q;„ - A)){{xl - OLn)Pn{x) - P„+l(x)), 

and the first-order differential equation (|3.29p : 

(A - Q;„ - (A - 0Ln)0Ln - CXn{A - OLn))Pn{x) 

= Pn{x)J - JPn{x) - x{Pn{x)A^ - A^ Pn{x)) 

+ {xl - A + an + (A - an)an + an{A - q;„))(P„(x)A - AP„(x)) 
+ (2/3„ + n(2A - /) - 2{A(3n + /3„A))P„(rc). 

Multiplying this equation by 2 and plugging it into the second-order differential equation ()4.8p 
gives 

P"(x) + 2P^(x)(A - xl) + P„(x)(A2 - 2xA^ - 2J) = (A^ - 22; A^ - 2j)P„(x) 
+ (2n(2A - I) - 4(A/3„ + /3„A) + 2((a„ - A)a„ + - A))A)P„(x) 

- 2((a„ - A)q„ + - A))(P^(x) + P„(2;)A). 

This differential equation is not of Sturm-Liouville type considered by Duran and Griinbaum 
in [23], but nevertheless we have been able to give a number of differential equations of first and 
second order satisfied by MOPRL with respect to a weight matrix that have not been considered 
up to this point. 

4.1.2 U'{x)U-^{x) = 2Bx 

The weight matrix ()4.ip is now given by 

where we assume that B is chosen such that all the moments exist. By (j4.2p . 
T{x) = e-^'/2gBx2 ^j^^ ^ - I)x. 

The weight matrix is an even matrix function, so that all moments of odd order vanish. As 
a consequence, a„^„_i = and Q!„ = 0. Then 

An{x;G) = 2{I - B* - -fnB-f~^), Bn{x;G) = {I-2B)x, 
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and there will be only one compatibility condition: 

2{I-B- 7-^iS*7„+i)/3„+i -2(3n{l-B- %^,B*-fn^i) = I. 
From (j4.16p we get the explicit expression for /3„ via 

2{I-B- 7-iB*7„)/3„ = nl + 2{an,n-2B - Ba„,„_2). 
The lowering and raising operators are 

P^(2;) + 2x{Pn{x)B - BPn{x)) = 2{I-B- 7-iS*7„)/3„P„„i(x), (4.16) 

and 

P^(x) + 2x{Pnix)B - BPn{x)) =2{I-B- 7^'-B*7„) (xP„(x) - P„+i(x)) , 

respectively. This yields the following second-order differential equation 

P^'(x) + 2xP^(x)(2S - /) + 2x(7-^S*7„ - i„)P^(x) + Ax^Pn{x)B^ - B^Pn{x)) 
+ AKnPn{x) + ((2 - 4x2)/ + 4x2(7-ip*7„ - L^)) {Pn{x)B - BPn{x)) = 0, 

where 

Ln = {l-B- -f-'B*-fn)B{l -B- j;^'B*-fny\ 

and 

Kn={l-B- -f~'B*Jn)f3n{l " P " 7„"^S*7n~l). 

These expressions are valid for any B (as long as all moments of the weight matrix exist). For 
further simplifications we can assume again that the hypotheses of Proposition 13.31 hold, and 
consider the cases given by the algebraic relations (|4.12p and (|4.15p . For instance, when ()4.15p 
holds, we obtain again a second-order differential equation with coefficients in the right hand 
side independent on n, studied by Duran and Griinbaum in [23]. The computations are similar 
and will be omitted for the sake of brevity. 

4.1.3 Other cases 

Assume now that we have a linear combination of the previous two cases, i.e. 

U'{x)U~^{x) = A + 2Bx, A,BeC^''^, with (7(0) = /, (4.17) 

in which case the matrix G{x) is given by G{x) = A + [2B — I)x. This is all that is needed to 
calculate the coefficients An{x; G) and Bn{x\ G) and consequently the compatibility conditions, 
the ladder operators and the differential relations. 

However, for the corresponding orthogonality weight we need to solve ()4.17p explicitly, which 
may be non-trivial (unless A and B commute). In general, this solution can be given in terms 
of the time ordered exponential 

U{x) = : er(^+s-)<^- :, A, P G C^>^^, x G M. 

This is obtained by rewriting the differential equation as an integral one and "solving" it by 
iteration. In general, one cannot give an explicit expression for this infinite sum. 

An example of explicitly solvable non-trivial equation (|4.17p can be found in j22|, Theorem 1.1], 
where non-commuting A and B are constructed (as a linear combination of the matrices L and J 
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in (|4.10p and ()4.1ip . respectively) such that the corresponding U{x) = Qi^-'"oJ)x^voJx ^ with vq 
any real number. [22] contains also another example of a pair of non-commuting matrices A 
and B such that U{x) is a matrix polynomial of degree A^ — 1, whose coefficients can be obtained 
recursively in terms of A and B (see Theorem A.l therein). 

An alternative way of generating examples is by considering weight matrices W of the 
form dm with U{x) G such that U{0) = I and det U{x) = 1 for x G M. If we denote 

m m 

C/(x) = / + J] AfcX^ U-\x) = I + ^Bkx\ 

k=l k=l 

then the coefficients Ak and B^ are connected by simple algebraic relations. For instance, 

Bi = -Ai, B2 = -A2 + Al, B3 = -A3 + A2A1 + A1A2- Al. 
In consequence, 

U'ix)U"\x) = Ai + (2A2 - Al)x + (3A3 - 2A2A1 - A1A2 + Afjx'^ 

+ (4A4 - 3A3A1 - A1A3 + 2A2AI + A?A2 + A1A2A1 -2Al-Af)x^ + ---. 

Let us discuss some cases, depending on the degree m of the matrix polynomial U{x). In the 
first non-trivial example, when m = 1, we have U{x) = I + Aix and U~^{x) = I — Aix with 
A1 = 0. Therefore, U' {x)U~^{x) = Ai, which was considered already in Section [4.1.11 

Let m = 2, so that U{x) = I + Aix + A2x'^ and U~^{x) = I - Aix + {A\ - A2)x'^. The 
condition U{x)U~^{x) = I yields the following algebraic relations between the coefficients Ai 
and A2: 

A2A1 + A1A2 = Af, A2(-A2 + A?) = 0. 

Algebraic manipulations of these two equations give additionally A2A1A2 = and A2AI = 
A\A2. Therefore 

G{x) = Ai + (2A2 - Al- I)x - A2Aix^. 

The simplest case when G{x) is a matrix polynomial of degree one gives another algebraic 
equation, A2A1 = 0. This immediately implies that A| = and A\ = 0, i.e. once again 
resulting in very strong algebraic conditions on the coefficients. In order to find an interesting 
example (when adyij(A2) 7^ 0) we have to go at least to the dimension = 4; for instance, 

/O 614 \ 
023034 
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Considering higher degree examples will give more and more algebraic relations among the 
coefficients A^, which in general are difficult to solve. 

In particular, in the examples that have appeared in the literature so far, U{x) is in general 
a matrix polynomial of degree depending on the size N . 
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4.2 The Freud case 

4 

Now we have q{x) = so the weight matrix is given by 

W{x) = e-''^U{x)U*{x), xGM. 
Again, we will consider two cases, first U'{x)U~^{x) = A and then U'{x)U~^{x) = 2Bx. 

4.2.1 U'{x)U-^{x) = A 

The weight matrix (j4.1|) in this case is given by 

W{x) = e-^■'e^^e^*^ A G C^^^, x G M. 
Then T{x) = e~^'*/^e"^^, G is given by 

G{x) = -2x^1 + A. 
Therefore using (j3.15p and (j3.16|) we get 

Anix; G) = 4{x'^I + OinX + /3„+i + /3„ + a^)*, 

and 

Bn{x\ G) = 2x^1 + 4(/3„x + I3nan-i + anPn) - A. 

These formulas are obtained by using the relation between the coefficients of the three-term 
recurrence equation ()2.30p and the coefficients of the MOPRL. In this case the coefficients 
Ani-',G) and Bn{-',G) are matrix polynomials of degree 2 and 3, respectively. 
From Proposition 13.51 the compatibility conditions are 

I + Aotn - OLnA - 4:{f3n+lOLn + CXn+lf^^n+ljOin + 4Q;n (/3„Q!„„i + Q;„/3.„) 

= 4(/3„+2 + /3„+i + q2_^i)/3„+i - AfinifBn + /3n-l + a^i), 

and 

4((/3„+i + (3n + a-DcXn + OCniPn+l + At) + ^n+lAt+l + PnOLn-l) = A + 'y~^A*Jn- 

The lowering operator is 

P^(2;) + Pn{x)A - APn{x) = -4(/3„X + PnCtn-l + CXnf3n)Pn{x) 
+ 4{x'^I + CXnX + I3n+1 + At + Q^) /9n-Pn-l (x) , 

while the raising operator is 

P^{x)+Pn{x)A- APn{x) 

= 4(x^I + (3n+lX - {Pn+1 + (3n + OlDan - f^nO-n-l " CXnfin) Pn{x) 

)Pn+l{x). 

From the term in the lowering operator we obtain 

An expression of the second-order differential equation can be given using Proposition 13. H 
but in this case it is necessary to compute the inverse of An{x\ G), which is a matrix polynomial 
of degree 2 with a rational function determinant in general. 
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4.2.2 U'{x)U-^{x) = 2Bx 

The weight matrix (j4.ip is given by 

we have T{x) = e~^*^^e^^' and 
G{x) = -2x^1 + 2Bx. 

Again, this weight matrix is an even matrix function, and we conclude as above that an,n-i = 
and cXn = 0. The coefficients (j3.15p and (j3.16p are given by (using the same approach as before) 

An{x; G) = A{x^I + /3: + /3:+i) - 2{B* + Tn^Tn ') , 

and 

Bn{x;G) = 2{x^I + {2l3n-B)x). 
Again, there is only one compatibility condition ()3.23p : 

I = 4((/3„+i + /3„+2)/3n+l - /3n(/3n + Pn-l)) 

- 2(S/3„+i - pnB + 7~+lS*7n - 7n '^*7n-l) 

with the lowering and raising operators 

Pl,{x) + 2x{Pn{x)B - BPn{x)) + 4x/3„P„(x) 

= (4(x2j + (3n+l +(3n)-2{B + -f-^B*-fn))f3nPn-l{x), 

and 

P^(x) + 2x{Pnix)B - BPn{x)) = (4x3/ + 2(2/3„+i -B- B* -fn) x) Pnix) 
+ (-4(x2j + /3,+i + /3„) + 2(S + 7-iS*7„))P„+i(x), 

respectively. 

Finally, considering the 0{x^~'^) term in the lowering operator we obtain 

nl + 2{an,n-2B - Ban,n-2) = 4(/3„/3„_i + f3l + /3„+i/3„) - 2{B + j~^B*-fn)Pn- (4.18) 

This equation may be regarded of as a matrix version of a discrete Painleve equation. Observe 
that in the scalar situation {B = 0), the equation (j4.18p reduces to the well-known string 
equations 

n = 4/3„(/3„+i+/3„ + /3„_i), 



studied in the context of discrete Painleve equations and orthogonal polynomials of Freud type 
in [52] (see formula (19) therein). 
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5 Final remarks 

This paper extends to the matrix case the methodology of derivation of differential relations 
for MOPRL using a RH approach. We obtain explicit formulas for the ladder operators of very 
general weight matrices factorized in the form W = TT*. We have also shown that in the 
matrix case there is an extra freedom absent in the scalar situation, which allows us to obtain 
a family of ladder operators, some of them of the 0-th order, which does not occur in the scalar 
case. Furthermore, by combining appropriately the family of ladder operators we get a family 
of second order differential equations satisfied by the MOPRL, some of them of the 0-th or first 
order (under some invertibility assumptions) . This yields some new results even in the particular 
cases studied before in the literature. 

In order to keep the size of this paper reasonable, we have restricted our attention to the 
weights supported on the whole real line. When supp(W) = [0, +oo) or supp(VF) = [—1,1], we 
can follow the ideas exposed in Section [3Tl except that now we have to assume that either zG{z) 
(for the case of the semi-axis) or (1 — z^)G{z) (for the case of the finite interval) are matrix 
polynomials. For the discussion of these situations in the scalar case, where the differential 
equations for the Laguerre and Jacobi polynomials are obtained, the interested reader is referred 
to [lOl Chapter 22]. 
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